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ABSTRACT 



Aims. We study the molecular content and chemistry of a circumstellar disk surrounding the Herbig Ae star AB Aur at (sub-)millimeter 
wavelengths. Our aim is to reconstruct the chemical history and composition of the AB Aur disk and to compare it with disks around low-mass, 
cooler T Tauri stars. 

Methods. We observe the AB Aur disk with the IRAM Plateau de Bure Interferometer in the C- and D- configurations in rotational lines of 
CS, HCN, C2H, CH3OH, HCO^, and CO isotopes. Using an iterative minimization technique, observed columns densities and abundances are 
derived. These values are further compared with results of an advanced chemical model that is based on a steady- state flared disk structure with 
a vertical temperature gradient, and gas-grain chemical network with surface reactions. 

Results. We firmly detect HCO+ in the 1-0 transition, tentatively detect HCN, and do not detect CS, C2H, and CH3OH. The observed HCO+ 
and ^^CO column densities as well as the upper limits to the column densities of HCN, CS, C2H, and CH3OH are in good agreement with 
modeling results and those from previous studies. 

Conclusions. The AB Aur disk possesses more CO, but is less abundant in other molecular species compared to the DM Tau disk. This is 
primarily caused by intense UV irradiation from the central Herbig AO star, which results in a hotter disk where CO freeze out does not occur 
and thus surface formation of complex CO-bearing molecules might be inhibited. 

Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks - individual: AB Aur - Radio-lines: stars 



1. Introduction 

The rich variety of the detected exoplanetary systems cannot 
be fully understood without knowledge about the physical and 
chemical evolution of their precursors - protoplanetary disks. 
The question what role chemistry plays during planet forma- 
tion and how it is coupled to disk dynamics remains loosly 
constrained. 

Nowadays it is widely believed that disk evolution is con- 
trolled by redistribution of angular momentum due to tur- 
bulent viscosity. A promising source of turbulence is the 
magnetorotational inst ability that works if the disk matter 
is sufficiently ionized ('Balbus & Hawley, '1991*). The chem- 
istry of ionization fraction in protoplanetary disks was stud- 



ied in detail, both theoretically and observationally (see 



e.g., G 


ammiel 19961; Glass sold et al. 


1999; 


Semenov, Wiebe. & Hennins 


20061: 


Pascucci et al.L 2007b. To a 



20041: lllgner& Nelson . 



nated by stellar radiation in the upper disk region and high- 
energy cosmic ray particles in their midplanes. For the inner, 
planet-forming disk regions, X-ray radiation from _a_young star 
may play a crucial role. Recently, iTurner et"aD (12007) inves- 
tigated the evolution of the ionization degree in the inner re- 
gion of a protoplanetary disk, using a coupled 3D magneto- 
hydrodynamical and chemical code and found that dynamic 
transport increases electron concentration within the dead zone. 
It has been widely discussed that the interplay between chem- 
ical and transport processes must have been important during 
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the formation of our own solar system (e.g. , Cvr et al 
Bockelee-Morvan et al. [ |2002l:IWooden et al.ll2QQ5k 



1998 



Due to high computational demands, the studies of disk 
chemistry coupled to disk dynamics are in their early stages 
Ilgner et al. [ l2004 IWillacv et al.L I2OO6I: iSemenov et al.L I2Q06I: 
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Table 1. List of observations using the Plateau de Bure Interferometer. In the column of the configuration, the abbreviation of 
e.g. "5D" refers to the number of antennas and the used configuration. 



Line transition frequency configu- baseline integration synthezied beam 1 a rms 

[MHz] ration range [m] time (Ton) [h] size, position angle [mJy/beam] 



HCO+ 




1-0 


89188.523 


5D,C2 


24- 


175 


5.9+10 


5.2'' X 4.8' 


,6° 


18 


HCN 


J=l- 


-0, F=2-l 


88631.847 


6CD 


24- 


175 


5.9 


6.0" 


X 4.2", 


105° 


24 


CS 




2-1 


97980.950 


5D 


24- 


160 


1.8 


9.0" 


X 4.2", 


120° 


170 


C2H 


N = 1-0, J 


= 3/2-1/2, F=2-l 


87316.925 


6D 


24- 


113 


1 


9.8" 


X 3.9", 


110° 


37 


CH3OH 


15(3,13) 


-14(4,10) A+ 


88594.960 


6CD 


24- 


175 


3.9 


6.1" 


X 4.2", 


105° 


24 


CH3OH 


15(3,12) 


-14(4,11) A- 


88940.090 


5D 


24- 


175 


3.9+2.2 


6.3" 


X 4.3", 


105° 


24 


CH3OH 


2(1,1) 


-1(1,0)^" 


97582.830 


5D 


32- 


73 


1.8 


9.4" 


X 4.3", 


148° 


170 



IXscharnuter & Gail 12007 ). Chemical models were success- 
fully applied t o study disk ch e mistry for steady- state accretion 
disks (e.g., Iwillacv et all 19981: ^an Dishoeck & Blake, 



cal models. In this second paper, we investigate the molecu- 



19981: lAikawa & Herbst . Il999[ Markwick et al.L 120021: 



Semenov et al .Dutrey et al.ll2007h . The most important 

result is the chemical stratification in such disks, where 
many molecules are abundant at an intermediate, slightly UV 
irradiated layer, while being frozen out in the cold and dark 
midplane and broken apar t by unshie l ded high-energy radia- 
tion in disk atmospheres. Bergin et al. ( 2003) have shown that 
the non-thermal UV radiation from T Tauri stars may partly 
come in Lyc^ photons, which significantly affects abundances 
of molecules dissociated by these photons. In contrast, Herbig 
Be/Ae stars emit much stronger thermal UV radiation, with 
circumstellar disks being hotter and more ionized than the 
disks around T Tauri stars. 

A detailed understanding of physical and chemical con- 
ditions of the planet-forming environment requires not only 
sophisticated models, but also high-resolution observations 
of protoplanetary disks, both in dust continuum and various 
molecular lines. Unfortunately, such studies are challenging 
due to the limited sensitivity and spatial resolution of avail- 
able obser vational facili t ies, a nd s mall disk size s 1 00 — 
1 000 AU). butrev et al.l d 1997b and lKastner etaP (Il997h have 
first detected several molecular species towards protoplane- 
tary disks by using the I RAM 30m single-di s h tele s cope, fol- 
lowed by observations of van Zadelhoff et al.l ( 20Ql ): lThi et al.l 
(l2004: iSemenov et al ]ll2005h. In terferometric observation s by 
[oil (l2QQlk : lAikawa et al.l (l2QQ3k : IPietu et al.l (120061 12007 ) al- 
lowed to study disk gas at outer radii > 50 AU in lines 
of several abundant species (CO, ^^CO, C^^O, HCO+, CS, 
CN, H2CO). These studies have confirmed that photochemistry 
plays an important role even at large disk radii and that abun- 
dances of many gas-phase molecules are depleted. 

In 2005 a joint Heidelberg-Bordeaux "Chemistry In Disks" 
(CID) project was established to investigate the spatial distri- 
bution of various molecular tracers across well-studied T Tauri 
and Herbig Ae disks of various age, followed by compre- 
hensive modeling. In the first paper of this series, we have 
presented the results of a deep search for N2H+ and HCO+ 
towards two T Tauri stars (DM Tau, LkCa 1 5) an d one 
Herbig Ae star (MWC 480), see butrev et al.l (l2Q07b . The 
N2H+ emission has been detected in LkCa 15 and DM Tau, 
with the N2H+ to HCO+ ratio similar to that of cold dense 
cores and the disk ionization degree as predicted by chemi- 



AB Aurigae ( 


van den Ancker et al.L 1998: Gradv et al.. 


Roberge et al. 


. 2OOII: iFukasawa et al.L l2004 distance 



pc), which is an intermediate-mass analog of the young T Tauri 
stars and evolutionary precursor of the main- sequence debris 
disk sources like P Pic, a Lyr, and a PsA. This disk is peculiar 
in comparison with other known circumstellar disks since it is 
still embedded in a rath er extended enve l ope a n d has clump\ 



densi t y sub- structures dFukagawa et al.L I2QQ4I: ICorder et 



2QQ5l:IPietu et al.ll2QQ5l: Lin et al.ll2006|). This points to ward a 



2W 



low-mass companion or a giant planet (|Rodriguez et al 112007), 
a recent stellar encounte r, or imprin t of ear lier non- steady evo- 
lution (see discussion i mPietu et al ](l2005h. hereafter Paper II). 



Using the IRAM 30-m telescope, ISemenov et al.l (I2005L here 



after Paper I) have studied the molecular content towards this 
source at low resolution and derived basic parameters of the en- 
velope. They have derived a low disk inclination of about 17° 
(face-on or ientation) and low disk mass of about 0.013 Mq. 
Pietu et al.l (.2005) have resolved the disk structure in different 
CO isotopes using the Plateau de Bure Interferometer (PdBI). 
They have found that the disk is inclined by about 30° and its 
rotation departs from the Keplerian law, with an exponent for 
the rotation velocity as low as 0.41. The outer radius of the 
gaseous disk is close to ~ 1000 AU, its inner hole is at a radius 
of 100 AU in the dust emission, and at ^ 70 AU in CO lines, 
and the disk mass is ^ 0.02 Mq. 

The main aim of the present study is to search for C2H, CS, 
HCN, HCO+, CH3OH, and CO isotopes in the disk of AB Aur 
using compact interferometer configurations of the PdBI, to de- 
rive their column densities, and to compare these values with 
those from a robust chemical model and the disk of a cooler 
T Tauri star DM Tau. 

2. Observations 

Observations with the PdBI were carried out between February 
and August 2002 using the compact C and D configurations 
and single side band tuning. Table [T] lists the observational pa- 
rameters for the individual observing runs. 

The frequencies are t aken from the Colo gne Database for 
Molecular Spectroscopy dMtiller etal.ll200l[ 2005). The phase 
reference center is RA(2000) = 04^55^54.8^ DEC(2000) = 
+30° 33' 04.3'.' All line observations were performed in the 
lower side band with a frequency resolution of 39 kHz, lead- 



Schreyer et al.: Molecules in the AB Aur disk 



30°33'10' 



30°33'00' 



^ 30°33'10' 

o 

o 

o 

^ 30°33'00' 



30°33'10' 



30°33'00' 



total 




5.09 




6.34 




4"55'"46^ 



4.09 



5.34 




6.59 




4.34 



5.59 





4.59 




5.84 




7.09 




4"55'"46^ 



4''55"'46^ 
RA(2000) 



4"55'"46^ 



4.84 



• 



6.09 




7.34 



(2) 



4"55'"46^ 



Fig. 1. Mapping results of the combined HCO+ (1-0) data. Contour levels are: -0.04 to 0.13 Jy km s~^ in step of 0.04 Jy km s~^. 
The mean velocity of each channel is given in each box top left. In the panel top left, the total integrated HCO+ map is shown. 
Contour levels are: -0.015 to 7 0.10 Jy kms~^ in step of 0.015 Jy kms~^. The area of the HPBW of the synthezied beam is 
indicated as a grey ellipse in the panel bottom right. 



ing to a velocity resolution of 0.14 kms~^. We took advan- 
tage of the flexible correlator to observe simultaneously transi- 
tions of CH3OH in the same bandpass as for the observations 
of HCO+, CS, and HCN. The total integration time was rather 
short, about 2 hours for each line. We combined the HCO+ 
(1-0) data obtained during these observations with that from 
Paper II to reconstruct the final image at 5.34'' x 4.8'' spatial 
resolution. 

The bandpass and phase calibration were performed using 
observations of the objects 3C84, 3C345, MWC349, CRL 618, 
and 0528+134. The GILDAS software is utiHzed for the data 
reduction and final phase calibration. The resulting synthesized 
half-power beam widths are listed in Table [T] 

All listed lines in Table [T] were centered to the systemic ve- 
locity of Vlsr = 5.8 km s~^. The weak underlying continuum 
was detected, with an integrated A 3 mm flux of 3.7 mJy. Thi s 
is consistent with previous measurements of Piet u et al.l (l2005b . 
For the final analysis of the observed line spectra this contin- 
uum was subtracted. 

3. Results 

The transition of HCO+ (1-0) was the only spectral line that 
was well detected in each observing run. In Figs. [T] and [2] we 
show the combined spectrum of HCO+ (1-0). With the beam 
of 5.34" X 4.8'' that is comparable to the disk size, we do not 
fully resolve the disk around AB Aur (Fig.[Tl upper left panel). 



However, the pattern typical for rotating circumstellar disks ap- 
pears clearly in various channel maps (Fig.[T]). The peak inten- 
sity of the double line profile of HCO+ is 0.12 Jy/beam, which 
is 10 times higher than the noise level. Fig. [2] shows the varia- 
tion of the line profiles across the spatial extent of the AB Aur 
disk, which is similar to the spectra presented in Papers I/II. 
The symmetric double-peaked line profile at the center of the 
spectral map is a typical feature of barely resolved Keplerian 
disks. Since deconvolution is a highly non-linear procedure that 
is prone to errors we perform our analysis in the iiv-space, 
using the -minimization technique of .Guilloteau & Dutreyl 
(Il998h . 

None of the other attempted species were detected (C2H, 
CS, and CH3OH), but we have a hint of the HCN (1-0) hyper- 
fine line components F = 1-1 and F = 2-1 at the 2 a noise level. 
The 1 cr rms noise levels for all species are listed in Table [T] 

To better quantify the column densities and molecular 
abundances, we analyze the data by applying the parametric 
disk model and the -min imization method as described in 
iGuifloteau & DutrevI (119981) . First, we analyze the HCO+ data 
leaving as many free parameters as possible. This allows us to 
verify the compatibility of the HCO+ emission with the disk 
parameters derived from the CO isotopologues in Paper II. The 
disk parameters derived from Paper II and found from HCO+ 
are presented in Tabled 

Although the angular resolution is low, the HCO+ data con- 
firms one important point: the existence of the inner hole of 
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:^ 75 — 100 AU radius found in Paper 11. This hole manifests 
itself as a lack of emission at projected velocities beyond ±1.6 
kms~^ from the systemic velocity, but the hole is not directly 
visible in the spectral map. The temperature is not well con- 
strained from HCO+, but compatible with that found from CO 
isotopologues. Overall, the non-Keplerian motion detected in 
Paper II seems less obvious in the HCO+ data: the best fit in- 
clination is slightly lower, giving a rotation velocity more in 
agreement with the stellar mass. Since our low angular reso- 
lution data masks out HCO+ emission from inner, highly in- 
homogeneous disk regions, the corresponding velocity profile 
seems to be closer to a Keplerian law. However, if the inner ra- 
dius is fixed to 75 AU, the velocity index derived from HCO+ 
becomes 0.36 ± 0.04, in agreement with Paper II. 

We use derived disk parameters (Col. 2 in Table O to deter- 
mine column densities for all observed molecules. The derived 
column density of HCO+ and upper limits for other species are 
summarized in Table (H These column densities scale linearly 
with the assumed temperature, but are independent from any of 
the other parameters. 

4. Modeling 

As a next step, we try to reproduce the measured column 
densities and the upper limits altogether using an advanced 
disk physical model and the chemical structure from a robust 
chemical model, similar to Paper I. The overall iterative fit- 
ting is easier than in Paper I since we refrain from comparing 
deconvolved observed and simulated spectral maps and base 
our analysis on the observed quantities derived with the x^- 
minimization method in the wv-plane. This allows us to ex- 
clude computationally expensive radiative transfer modeling 
with beam convolution from consideration. 

To simulate the disk phys ical structure we utilize a pas- 
sive flared 2D disk model of DuUemond & DominikI (2004) 
with a vertical temperature gradient and non-grey dust opac- 
ities. Some input disk and stellar parameters are taken as 
determined by the -fitting or found in pre v ious works 



van den Ancker et al.L . 1 998.: iThi et al.l 120041: IPietu et al 



20051: ISemenov et al.l Eoosb . see Table [H Other key parame 



ters, like the radial slope of surface density, the disk mass, and 
the cosmic ray ionization rate, are found by iterative fitting of 
the observed data (Table O. 

We assume that the central AOe star has an effective tem- 
perature ofj^^00^2jTadiusof2^^^ a mass of 2.4 solar 
masses (I van den Ancker et al. The dust grains are 

modeled as compact spheres of uniform 0.12 /im radius made 
of pure amorphous s ilicates with the optical data taken from 
Draine & Led (Il984|) . The standard 1% dust-to-gas mass ratio 
is used. Utilizing the results of Papers I and II, the disk inner 
radius is assumed to be ~ 70 AU, the outer radius is 1 100 AU, 
and the disk age is between 2 and 5 Myr. 

The disk is illuminated by the UV radiation from the cen- 
tral star and by the interstellar UV radiation. The inte nsity of 
the un -attenuated stellar UV flux is calculated using the lKurucz 
(19931) ATLAS 9 of stellar spectra and converted to the x* = 
lO^X factor at 100 _AU^wherethe factor x is the mean inter- 
stellar UV field of lPraind (|l978|) . The UV intensity at a given 



Table 2. List of disk and stellar parameters using the x^ mini- 
mization in the UV plane 



disk parameters derived from 
CO (Paper II) HCO+ (this paper) 



Vlsr (kms-^) 


5.87 


5.87 ±0.03 


Inclination (°) 


35 


21 ±16 


Vsin(i) (kms"^) 


1.67 


1.49 ±0.07 


Velocity exponent 


0.40 


0.40 ±0.04 


Position angle (°) 


-30 


-20 ±4 


Inner radius (AU) 


75 


110 ±16 


Outer radius (AU) 


1000 


> 800 


Temperature (K) 


35 


[35] 


exponent q 


0.1 


[0.1] 


Column density (cm~^) 


6.2 X 10^^ 


6.2 ±1.2 10^^ 


exponent 


2.5 


3.5 ±0.5 


Line width (km s ~ ^ ) 


0.40 


0.23 ±0.10 



The fixed parameters are taken from Table I of Paper II, as an average 
of the results from the ^^CO. The column density is for H2 (Col 2) 
and HCO+ (Col 3) at 250 AU. 



Table 3. List of disk and stellar parameters used in the best-fit 
model 



Parameter 


Dimension 


Fixed 


Derived 


Stellar radius 


Rq 


2.5 




Stellar mass 


Mq 


2.4 




Effective temperature 


K 


10000 




Stellar UV RF 


Xo 


100000 




at 100 AU 








CRP ionization rate 


s-i 




410-^« 


Inner disk radius 


AU 


70 




Outer disk radius 


AU 


1 100 




Disk mass 


Mq 




0.02 


Surface density 






profile 






-2.15 


Grain size 


/im 


0.12 





Dust-to-gas 
mass ratio 



0.01 



The fixed stellar and disk pa ramet ers are taken from lPietu et al.l 
(l2005h . 



disk location is calculated as a sum of the stellar and interstel- 
lar components that are scaled down by the visual extinction in 
vertical direction and in direction to the central star (ID plane- 
parallel approximation). We m odel the attenuation of cos- 
mic ra ys (CRP) by Eq. (3) from Sem enov. Wiebe. & HenningI 
(I2OO4I) and vary the initial value of the ionization rate Ccrp be- 
tween about 10~^^ and 10~^^ s~^ to match the HCO+ data. 
In the disk interior, ionization due to the decay of short-living 
radionuclides is t aken into account, assuming an ionization rate 
of 6.5-10-^^ s-i (lFinocchi&GailLll997h . The stellar X-ray ra- 
diation is assumed to be weak in intermediate-mass stars due to 
a lack of the dynamo mechanism and thus is neglected. Though 
the observed X-ray luminosity of the AB Aur is non-negligible, 
its spectru m is dominated by soft photons with energies below 



1 keV (.Telleschi et al.L 120071 which cannot penetrate easily 



deep into the disk. 
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offset [arcsec] 




2 4 6 8 2 4 6 

Velocity [km s '] 
Fig. 2. Spatial variation of the HCO+ (1-0) line profiles across the AB Aur disk. 



The gas-grain time-dependent chemical model adopted in 
this study is mostly the same as in Paper I, but with two ma- 
jor modifications. First, we now use osu.2007, the latest ver- 
sion of the O hio State Univers ity (OSU) database of gas-phase 
reaction^ ( Smit h et all 12004 ). with all recent updates to the 
reaction rates. Second, we utilize a standard rate approach to 
the surfac e chemistry mode ling but without H and H2 quantum 



tunneling (Katz et al.Lll999l) . So, only thermal hopping of sur- 



face reactants is allowed. The surfa ce reactions together with 
desorption energies are taken from iGarrod & H erbst ( 2006), 
which is mostly based o n previous stu dies of Hasegawa et al. 
(|l992|) and Hasegawa & HerbstI (Il993l) . Finally, rates of sev- 
eral tens of photodiss ociation and photoionizatio n reactions are 
updated according to ^ van Dishoeck et aP (2006), using the UV 
spectral shape typical of a Herbig Ae star and the ISM-like dust 
grain optical properties. 

Overall, the disk chemical network consists of about 650 
species made of 13 elements and 7300 reactions. Using this 
time-dependent chemical model, distributions of the molecular 
abundances and column densities for the considered species are 
simulated over 5 Myr of evolution. 

5. Results and Discussion 

The Fig. [3] shows the calculated 2D abundance distribution 
for all observed molecular species. The vertical scale is given 



See: https://www.phy sic s. ohio-state. edu/ eric/research, html 



in units of p ressure scale height calculated for t he midplane 
temperature (Dartois, Dutrey. & GuilloteauLl2003h . The corre- 
sponding vertically integrated column densities versus the ra- 
dius are shown and also listed in Tabled This best-fit physico- 
chemical model of the AB Aur disk is obtained after about 20 
iterations by varying the cosmic ray ionization rate, total disk 
mass, and surface density exponent. All values refer to a radius 
of 250 AU. Upper limits for CS and CH3OH were derived using 
the 2 (J errors of the -minimization. For these estimates, the 
H2 column density of A^(H2) = 6.2 x 10^^ cm~^ was adopted 
from Paper II as a reference. We have adopted an ^^C/^^C iso- 
topic ratio of 70. 

In agreement with previous studies and the observations of 
T Tauri disks, chemical stratification is apparent in the disk 
around AB Aur (see Fig. 5 in Dutrey et al.. 2007). The disk sur- 
rounding the central AOe star is hotter and has a warm midplane 
(T > 20 K) and, therefore, CO molecules do not freeze out 
while other molecules do. Due to self- shielding, energetic UV 
radiation from AB Aur cannot dissociate a large fraction of CO. 
This leads to a high CO column density of about 3 10^^ cm~^ 
at 250 AU, with a relative CO abundance of 7 10~^. It seems 
that such a low surface CO population in the disk of AB Aur 
may lead to inhibited catalytic formation of complex molecules 
like methanol via hydrogen addition reactions on dust surfaces. 
In contrast, in the DM Tau disk a substantial fraction of CO 
resides on dust grains, so that the averaged relative CO abun- 
dance in the gas phase is 6 times lower, ~ 10~^. We found that, 
within adopted disk masses of 0.001 — O.lMgun, the radial pro- 
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file of the CO column density obtained by chemical modeling 
strictly follows the input surface density profile with a scaling 
factor. 

The chemically related molecular ion, HCO+, is directly 
produced from CO by ion-molecule reactions with and gets 
destroyed by dissociative recombination. In turn, formation of 
H3" is solely due to interactions of cosmic ray particles with 
molecular hydrogen (see e.g., lOkaL [2Q06|) . Thus, the HCO+ 
column density is sensitive to adopted cosmic ray ioniza- 
tion rate and, to some extent, temperature structure. Our best- 
fit model reproduces the observed HCO+ column density at 
250 AU, 6 10^2 cni-2, within a factor of 2.5 (iVmod(HCO+) ^ 
10^^ cm~^), which is co mparable to observati onal and inherent 



chemical uncertainties dVasyunin et al.U2QQ8h . 

In contrast to the DM Tau disk where HCO+ is mostly con- 
centrated in the intermediate layer due to freeze out of CO, 
HCO+ reaches the highest concentration in the midplane of 
the AB Aur disk. The HCO+ ion was found to be very sen- 
sitive to the adopted value of the CRP flux and radial density 
distribution. The acceptable fit to observations was only pos- 
sible when the CRP ionization rate was lowered by a factor 
of 3-5 compared to the standard value of 1.310"^'' s~^, the 
surface density profile is ^ —2.15, and the disk mass is be- 
tween 0.008 — 0.012Msun. With lower disk masses, the pre- 
dicted ^^CO column density was lower than the observed value 
(4 10^^ cm~^ at 250 AU), while with other density profiles it 
was impossible to get proper CS, C2H and HCN column den- 
sities. The model results dependent on the adopted inner radius 
- with radii smaller than about 50 AU and fixed total disk mass 
the modeled HCO+ and CO concentrations were lower than the 
observed values. As we used the newly measured inner radius 
of ~ 70 AU (see paper II), it is not surprising that the mod- 
eled distribution of HCO+ in the AB Aur disk differs from the 
results of Paper I (see Fig. 5 therein). 

The reason why the CRP ionization rate could be lower in 
AB Aur than a widely accepted value 10~^^ s~^) is many- 
fold. First of all, the standard ionization rates of He and H2 can 
be uncertain by a factor of 2 (IWakelam et al. ['2006). Second, if 
a self-generation of magnetic fields in protoplanetary disks by 
a dy namo-like mechanism takes place ( Dolginov & Stepinski 
Il994|) . then the cosmic ray particles will be scattered and will 
not reach disk midplanes - exactly the zone where a significant 
fraction of HCO+ exists in the AB Aur dis k. Third, the AB 
Aur disk has a complex and clumpy struct ure (iFukagawa et al.. 
I2OO4I: IPietu et al.l l2005l: iLin et al.l l2006h . and, consequently, 
efficiency of dissociative recombination of HCO+ will be dif- 
ferent in clumps compared to an inter-clump medium. Our ax- 
isymmetric disk model does not take this effect into account. 
Last, there may be some missed important reaction pathways 
to HCO+ or improperly derived reactions rates. 

Abundances and column densities of other observed 
molecules are also explained by this model. The chemically ac- 
tive C2H radical is only abundant in dilute, irradiated disk sur- 
face and quickly converted into heavier, more complex species 
in the denser region adjacent to the midplane (Fig.O. On the 
other hand, HCN and CS are more widely distributed across 
the AB Aur disk, though both these molecules are underabun- 
dant in the midplane where they stick to grains and surface pro- 



cessed. Note that the column densities of CO and thus HCO+ 
decrease with radius, whereas the column densities of C2H 
and HCN slightly increase with radius. The CS column den- 
sity stays nearly constant for large radii, r > 200 AU, and is 
found to be stable against iterative variations of the input pa- 
rameters. It is HCN that depends on the adopted value of the 
surface density exponent. As this molecule can be easily pho- 
todissociated by stellar UV photons, its abundance is sensitive 
to details of the disk vertical structure, in particular, shielding 
by dust. Presumably, all modeled column densities agree well 
with the observed values or upper limits due to the lack of sen- 
sitivity of the measurements. On the other hand, the derived ra- 
dial profiles are shallower than extracted from observations by 
the -minimization analysis. The surface density is a measure 
of the total disk mass, which is usually obtained from thermal 
dust emission and poorly known dust properties (gas-to-dust 
ratio, opacities), and thus can be rather uncertain (see discus- 
sion in Paper I). 

It is interesting to compare the molecular contents of the 
AB Aur disk and of another well- studied disk around the T 
Tauri star DM Tau. The DM Tau disk ha s a similar size of y 
1 000 AU with a hole of a few AU in size dCalvet et al.l.l2Q05h . 
is more evolved (~ 5 Myr), and more massive (~ 0.05 Mq) 
(see Table 2 in Dutrey et al.Ll2QQ7|) . In Table|4]the observed col- 
umn densities reseated to the ^^CO value at 250 AU are listed 
for the DM Tau disk as well (last column). These reseated val- 
ues are less sensitive to possible errors in the estimate of the 
disk mass and surface density. 

It is apparent that relative column densities for all observed 
species are lower in the disk of AB Aur compared to the values 
observed in the DM Tau disk. Given that the absolute column 
density of ^^CO at 250 AU is similar in both disks, it implies 
a poor molecular content of the AB Aur disk (except of CO 
and H2). However, the disk around DM Tau is about 4 times 
more massive than the AB Aur disk. The fact that both disks 
harbor similar amounts of CO despite significant difference in 
their masses is likely a direct manifestation of severe freeze out 
in the cold midplane of the DM Tau source and/or difference in 
grain properties. 

The lower mass of the hotter AB Aur disk is likely the rea- 
son for its molecule-deficient gas content. Since the total dust 
shielding is lower, while the UV luminosity of the Herbig AO 
star is orders of magnitude higher than that of DM Tau, the pho- 
todissociation of many molecules is more effective in the disk 
around AB Aur. The age difference between the two sources is 
a negligible factor for the results of chemical modeling, even 
though gas-grain models usually do not reach a steady- state 
within a few Myr. The same is true for surface chemistry, which 
is not much of importance for most of the observed species, 
apart from methanol. 

6. Conclusions and Summary 

The AB Aur environment is clearly dominated by the envelope, 
as shown by the detections made with the 30-m telescope and 
the lack of detections with the IRAM array in short (1-3 hours) 
integrations. The disk has an inner hole in the gas distribution 
with a radius of ^ 70 AU. Using a gas-grain chemical model 
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-minimization method 


Chemical model 


DM Tau 




N 


1 (7 




N 








[cm-2] 


error 




[cm-2] 






H2 


610^^ 


110^^ 


1.510^ 


510^^ 


1.3 10'' 


110' 


i3co(*3) 


410^^ 




1 


410^^ 


1 


1 


HCO+ 


610^^ 


310^^ 


1.510"^ 


1.510^^ 


410-^ 


210"^ 


HCN 


510^^ 


310^^ 


1.310"^ 


410^^ 


10"^ 


710-^ 


CS 


310^2 


310^2 


< 810-^ 


210^^ 


510-^ 


310-^ 


C2H 


210^^ 




<510-^ 




2.510-^ 


10"^ 


CH3OH 





710^^ 


< 210"^ 












Relative to the ^^CO column density at 250 AU obtained by the -minimization method, 
^^*^ Relative to the ^^CO column density at 250 AU obtained by the chemical modeling, 
see results reported by Pietu et al. (2005). 



with surface reactions coupled to a flared passive disk model, 
we reproduce all observed column densities and upper limits. 
Note that such a comparison between the advanced chemical 
model and the interferometric observations are hampered by 
the lack of knowledge on the H2 surface density. The surface 
density of the AB Aur disk is derived from the millimeter dust 
emission, a procedure prone to uncertainties. Part of the dis- 
crepancies between the model and the observations may also 
be due to inhomogeneous structure of the AB Aur disk and/or 
grain evolution in central disk regions that is not considered in 
the calculations. Modeled and observed column densities rela- 
tive to ^^CO are both lower for the AB Aur disk than the values 
measured in DM Tau. The absolute amounts of CO gas are sim- 
ilar in both disks. The poor molecular content of the AB Aur 
system can be explained by much more intense UV irradiation 
by the central AO star and less effective dust shielding in the 
low-mass disk compared to the more massive disk around the 
less luminous Ml star DM Tau. The indirectly inferred cos- 
mic ray ionization rate that is needed to fit the HCO+ data is a 
factor of 2-5 lower than the standard value, which can be ex- 
plained by scattering of CRPs in a magnetized AB Aur disk or 
a clumpy disk structure. 
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Fig. 3. (Left) Modeled abundance distributions in the disk of 
AB Aur at 2 Myr (relative to total amount of hydrogen nuclei) 
for H2, CO, HCO+, HCN, CS, and C2H. The vertical axis is 
angle in radians. (Right) The corresponding radial distributions 
of the column densities (cm~^). 
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